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Abstract 

In this paper we examine an alternative formulation of the gauge principle in which the emphasis 
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the basic philosophy on mechanical systems, and then with the help of De Donder-Weyl formalism 
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Nambu-Goto closed string and relativistic field theory. 
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I. INTRODUCTION 



The gauge principle (see [1] for an historical account) is a basic ingredient of modern 
theoretical physics, particularly in quantum field theory. It is not necessary to elaborate 
much on this undisputable fact. A quick presentation of its main idea is that by gauging 
a rigid symmetry one must pay a "price" : that of introducing a new field, the gauge field, 
which geometrically represents principal connection on a principal bundle. This "price" 
has turned out to be an unexpected bonus which has irrevocably changed the theoretical 
landscape in physics. 

In this paper we propose to revisit the gauge principle from the point of view of enforc- 
ing constants of motion as constraints. Wc should, however, forewarn that our subsequent 
considerations will be purely classical, so particularly ordering issues will be outside our 
scope. Similarly we will assume that Lagrangian/Hamiltonian systems are equivalent when 
they produce identical equations of motion (EOM). This "on mass-shell" (i.e.. the classical 
path) identification is clearly not satisfactory from a quantum point of view where also "off 
mass-shell" behavior non-trivially contributes into, say, transitional amplitudes. Grassmann 
variables will also not be considered, since that complication is a straightforward general- 
ization. 

It is well known that theories — derived form a variational principle — which exhibit 
gauge invariance must be described by constrained systems. With these two words we refer 

to the framework put forward by Roscnfeld [2, 3], Dirac [4, 5] and Bergmann [6-8], who, 
independently, laid the ground to deal with such systems. In particular Roscnfeld's contri- 
bution, which has been overlooked for a long time, has recently resurfaced thanks to the 
work of D. Salisbury and it is discussed in [2]. The constrained systems are characterized 
by Lagrangians whose Hessian matrix with respect to the velocities is singular, thus pre- 
venting the Legendre map (LM) from tangent bundle (i.e., positions and velocities space) to 
cotangent (i.e., space of positions and momenta, or phase space) from being invertible. It is 
precisely the singularity of the Hessian matrix which makes room for the possible presence 
of gauge freedom. Eventually, the picture obtained in phase space is that we have a (non- 
uniquelly defined) canonical Hamiltonian - /J, and a set of primary constraints Ca that are 
just the consequence of the non-invertibility of the LM. Thus the dynamics in phase space 
is given by the Dirac Hamiltonian, 

Hd := H + X'^Ca, (1) 

with A" a set of in principle arbitrary Lagrange multipliers, together with the requirement 
that motions must satisfy the primary constraints, 

Ca = 0. (2) 

Here we will not dwell in the details of the theory of constrained systems, for which we 
simply refer to the literature [9-12]. What we want to emphasize is that, given the structure 

of the dynamics in phase space, one could think of a process of gauging a regular theory 
by just starting with an ordinary Hamiltonian H and a set of functions Ca that are to be 
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enforced as constraints. Then we could define a new dynamics by tfie equations (1) and (2), 
which hopefully would describe a gauge theory. 

In general, this program is bound to fail because the constraints must have a certain 
degree of compatibility with the generator Ho of the dynamics. Geometrically, one needs 
the dynamical trajectories to be tangent to the surface defined by the constraints. In general, 
one expects this requirement to eventually end up with the appearance of new constraints as 
well as the determination of some of the Lagrange multipliers. But if H and are chosen 
too arbitrarily, the most likely outcome is that there will be no set of A"'s that keeps the 
dynamical trajectories tangent to the constraint's surface. 

But there is a nice exception, with plenty of interest: if we choose the would-be constraints 
as some of the constants of motion for H, then full compatibility is easy to achieve. This 
is the case we will explore. We consider a Hamiltonian for a regular theory - obtained 
from a Lagangian in tangent space through an invertible LM - and a set of constants of 
motion satisfying {Ca, H} — and closing a certain algebra {C^, C;,} = c^j,Cc with c^j, 
being structure constants. For simplicity's sake we restrict ourselves to constants of motion 
without explicit time dependent, i.e. to scleronomic constants of motion. We then declare 
that the new dynamics is governed by the Dirac Hamiltonian (1) under the condition that 
the constants of motion are enforced now as constraints (2). 

To check that we are indeed on the right track we must verify that with these conditions 
the theory defined by (1) supports gauge symmetries and that they act on the "matter" 
fields as they should, just generalizing the action of the former rigid symmetries. Once this 
check is done, we can explore the new gauge theory and its dynamical consequences, because 
the dynamics is expected to undergo important changes after the gauging of the rigid group 
of symmetries. Finally we can further modify the theory in a natural way by introducing 
gauge invariant kinetic terms for the Lagrange multipliers. The full-fledged gauge theory is 
then obtained, with the new non trivial interaction terms allowed by the gauge principle. 

Our paper is organized as follows: In Section II we formulate our basic strategy using 
the language of mechanics. Namely, we show how to construct a gauge invariant theory by 
promoting constants of motion to constraints. We also stress an intimate connection with the 
mathematical structure of non-abelian Yang-Mills theory [13]. In Section III we complete 
the theoretical setup. Examples in mechanics are given in Section IV and the relativistic 
field theory is dealt with in Section V, where the key role of the De Donder-Weyl formahsm 
is made manifest. We devote Section VI to the case of the closed bosonic string and use our 
approach to obtain world sheet general covariance. Finally, we conclude in Section VII with 
a brief summary of our results and outlook. 

II. THE NEW GAUGE THEORY 

We start by considering a Hamiltonian for a regular theory together with a set of scle- 
ronomic constants of motion Ca satisfying {Cq,, H} = and closing an algebra {C„, C^} = 
c^(,Cc. Now we will prove that when Ca are enforced now as constraints then this new 
theory is indeed a gauge theory. The simplest way to prove it is by defining the extended 
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Lagrangian (indices for vector components are normally suppressed) 

Le{q,p,q,p,X) = pq- H{q,p) - X^Caiq^p) , 



(3) 



and proving that it has gauge transformations. Note first that the EOM for (3) coincide 
with those derived from the Dirac Hamiltonian (1) and the constraints (2) (an advantage of 
the Lagrangian formulation is that all the dynamics is encoded in a single function). Note 
also that we have enlarged the configuration space with the multipliers A" as new variables. 
We will prove that indeed (3) has Noether gauge symmetries. Since the constants of motion 
are the Noether generators of the rigid symmetries, it is reasonable to expect that the 
generator of the would-be canonical gauge transformations can be written as G = e"'{t)Ca, 
with e" being a set of arbitrary time-dependent functions. We will prove now that indeed G 
generates gauge transformations. The corresponding variations can be written as 



S^q' = {q\G}, 5ePi = {Pi,G}, 



(4) 



and the variations of the multipliers will be determined below by the condition that, under 
the variations thus defined, the Lagrangian is quasi- invariant, i.e.. 



for some F linear in e and its derivatives. Indeed, 

Se{Pi q') 



(5) 



Pi64' + q'^ePi = -PiSeq' + q'SePi + -^{piSeq') 



dpi 



(6) 



and 



5,{H + A«C„) = e\t){H, C,} + CaSA" + A"e''(t){C„, = C„(5,A« + X''e\t)cl,C, . (7) 



a ,6/ 



Thus, the appropriate definition 



{Doe{t)r, 



(8) 



makes S.Lp to be 



Ca 



d 
di 



6«(t) 



dCa 

' dpi 



Ca 



(9) 



which proves that the variations (4) and (8) define a Noether gauge symmetry for Lg. In 
Eq.(8) we have introduced the covariant derivative [14] 



(Do): ■■= m - X'cl, 



(10) 
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which is nothing but the covariant derivative for the adjoint representation. AnalogoTisly, 
one can introduce the covariant derivative for the phase-space variables = {pi, . . . ,q^, . . .} 
as 

DoC dtC - A«r„e = dtC - A^e, CJ (11) 

Here is the representation of the symmetry generators Ca that acts on Indeed, the 
Jacobi identity for Poisson brackets ensures that [F^, F^] = — cJ^^Fg. Because both dt and 
{ ,Ca} fulfil the Leibniz rule one can extend the covariant derivative (11) to any function 
0(£) on phase space. 

Using the fact that our active variations commute with the time derivatives, i.e., 5^{dt(f)) — 
dt{Se(f)), it is easy to check that the covariance condition takes the form 

d.iDo 4>) = e"i^o({0, Co}) =: e"F„(Do 0) . (12) 

The first equality in (12) can be proved by considering: 1) that active variations commute 
with the time derivatives, and so Fa((9t0) = Ca}, and 2) that the action of Fa on 

the multipliers is the adjoint action: FaA^ = c\^y^. Notice that in the last equality in 
(12) the representation F^ of the symmetry generators acting on D^q was defined. This 
definition turns out to be an exact identity on mass-shell. In this regard it is interesting 
to reahze that the curvature F — [Dq, Dq] = 0, and so in the case of mechanics the usual 
gauge-invariant kinetic term Tr(F^) is trivially zero. Thus the multipliers A" cannot become 
dynamical variables. On the other hand, a subsequent elimination of the momenta — which 
are auxiliary variables (auxiliary variables are by definition variables that can be isolated 
by using their own EOM) for Lg — , as done in the next subsection, will assign the A"'s the 
status of auxiliary variables. 

Note that (8), (10) and (11) carry indeed all the flavor of the transformation of a gauge 
fleld in a non-abelian gauge theory. This is exactly the case, because what we have done 
is precisely the application of the gauge principle: to gauge a group of rigid symmetries. 
We remind that the rigid symmetries are generated by constants of motion while the gauge 
symmetries by the first class constraints [15]. Thus gauging a group of rigid symmetries is 
tantamount to enforce the generating constants of motion as constraints. In this respect A" 
play the role of a connection in a principal bundle over M. The fact that F = then indicates 
that this bundle is flat (not big surprise for a bundle with so simple base space). Let us, 
however, stress that our derivation would go through even if not very rigid symmetry is 
gauged. For instance, we could have limited ourselves only to gauging any subgroup of rigid 
symmetries. The analogy with non-abelian Yang-Mills theory is summarized in Table I. 

Note finally that the case of a soft algebra [10] is easily accommodated. We can relax the 
condition that the constants of motion Ca form a Lie algebra to that of a soft algebra, where 
there are no longer structure constants but structure functions, {Ca, Ch} = c'^fj{q,p)Cc, and 
we can also relax the constant of motion condition, {Ca, H} — 0, to {Ca, H} = a\{q,p)Cb. 
In this case, equation (8) changes to 

8X := e'^it) - X'e\t)cl - e\t)al (13) 
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TABLE I: Comparison between the gauge theory presented in Section II and the non-abehan Yang-Mills 
theory. The parallelism obtained allows to formally identify Dq o Dfj, and A o A^. 



Gauge theory from Section IP 


Non-abelian Yang-Mills theory'' 




5A-{x) = D^e<'{x) 


= mm 


5^{x) = iT{e)^{x) 




e'^ix) = df.e'^ix) + A';^{x)f^,, e-(x) 


Dom^dtm-^''^am 


D^^{x) = - lAfJa^ix) 






On-mass-shcU situation {dtCa = 0) 


Doe = dte - {A, e} 


Df, e = d^e + [A^,e] 


(5,(I?o0) =i?o({</',e}) =r(e)Do0 


d{D^^) = iT{e)D^^ 



"Here we accept notations: A = A°Ca, e = f°'Ca, £. = {pi,P2, • • • , ,<P , . . .} is a phase-space point, </) is an arbitrary function 
on a phase space and T{Ca) = Ta = { , Ca} = a;*-' 

''Here we ax;cept notations: = —iA'^ta, e = e"ta, * is an arbitrary field multiplet and T(ta) = Ta is an irreducible 
representation of the algebra of ta generators that is adapted to ^, e.g. for ^ in fundamental rep. of SU{N) then T(Ajj) = 
—lA'^Ta with Ta being the [N x N) hermitian matrices. Generators in self-adjoint rep. are defined as, (T(,)° = if^^- 



In field theory one can find more general cases [10], like that of an open algebra, where 
the algebra of the constants of motion only closes up to linear terms that are antisymmetric 
combinations of the equations of motion, or when there is functional dependence among 
the constants of motion. We believe that these cases can also be addressed, but since the 
ordinary case already requires a non-standard formalism (see Subsection V C) we leave them 
for further study 



III. INVERTING THE LEGENDRE MAP 
A. The Lagrangian Lx 

Consider now the equations of motion for L^, i.e., 

C„ = 0. (14c) 
It is interesting to observe that by introducing the symplectic matrix u) 

= ( ^1 1 ' ^^^^ 
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(with Uj^j^ — u!^^) the EOM (14) can be succinctly written as 

Cn. = 0. 



(16) 



We can now use the first set of Eq.(14) to isolate the momenta in terms of positions q, 
velocities q and the multipliers A", thus rewriting (14a) in the equivalent form 



p - P{q,q,X) = 0, 



(17) 



for some functions P. This invertibility of the LM will hold in general. In fact, since the 
starting theory was not gauge, invertibility is guaranteed for A" — 0. With A" being just 
new independent variables, invertibility will be maintained in general. 
We implement p — > P(g, q, A) into to define the new Lagrangian Lx, 



Lx{q,q,X) = P{q,q,m - H{q, P{q,q, X)) - X^^C^iq, P{q,q, X)) ■ 
Notice then 



(18) 



dq 



V dq 



dq 

adCa 



+ 



V dq dq 



\ dp dv y 



adCa 



dp 



dP 

p^p 9q 



(19) 



(J - dH - \adCa 
" dp dp 



p^P 



because ^\p-^p = (q 

define the functions P{q, q, A). By the same token, we obtain 



vanishes identically owing to the procedure to 



dq 



^(?,?,A), 



(20) 



so we reobtain the functions P as the definition of the new Lagrangian momenta. By taking 
into account Eq.(20) the reader may note that the EOM for the Lagrangian Lx yield 



dP dH ^^dCa _ Q 
di dq dq 



(21) 



which is equivalent to (14b) when the identity (17) is utilized. 

The remaining EOM for is the one associated with the multiplier A. This equation 

sets just the constraint as EOM. From the perspective of Lx-, we can write down EOM for 

vanishes identically, we 



\a. OLx 

end up with 



dp Ip-s-P dX 



^ - Ca{q,P{q,q,X)), but again, since 

Ca{q,P{q,q,X)) = 0, 



dp \p-i-P 



(22) 



as the last EOM for Lx- This shows the equivalence between EOM from Lc and Lx, because 
(22) is just (14c) with the substitution p P{q, q, A), which is nothing but the EOM (14a). 



7 



B. Gauge symmetry for Lx 



Let us now prove that Lx has the gauge symmetry 6xq = {6eq)\p^p, 5aA = ^^A. One has 



SxLx = 



. dq 



Sxq + 



dq 



Sxq + 



PL, 
dp 



hP + 



p-^p 



5xX.(23) 



p-^p 



We need not care of the term with 5xP because (^)|p_j.p = identically due to the equiv- 
alence between (14) and (17). As regards Sxq we can write it as {Seq)\p-^p- All in all we can 
write 



5xLx = (5,L,)\p-,P = (^^F^ 



p^p 



dt 



(P \p^p) ' 



(24) 



where [17] in the last step we use that p ^ P implies also p — >■ ^P, etc. Thus we have 
proved that Lx inherits the gauge invariance of L^. 



C. A step further 

Finally, if A can be isolated from the equation (22), this means that it is in fact an 
auxiliary variable. It is well known that auxiliary variables can be substituted back into 
the Lagrangian without affecting the dynamics (see e.g., the Appendix in [16]). In fact 
the earlier substitution p — > P(gf, q, A) in the previous subsection is an example of this 
mechanism, for the variables p are isolated by use of their own equations of motion (14), but 
we have been explicit in the proof of equivalence of EOM. Thus with the substitution we 
would have arrived at a new Lagrangian L{q,q) with a dynamics equivalent to that of Lg. 
Of course, there may be technical obstacles to carrying out this step: solving the system of 
equations (22) may prove too difficult, getting rid of the multipliers can lead in general to 
impractically complicated, non polynomial expressions for L, etc. One can then revert back 
to the Lagrangian Le, with its EOM (14). 



D. • • • and a step further 



Despite potential complications related to solving the system (22) wc will suppose that 
indeed the variables A" can be isolated from the equations (22) and eliminated by plugging 
them back into the Lagrangian Lx- Thus (22) will be equivalent to A" = A"(g, q) for some 
functions A". We will prove that L{q, q) := {Lx{q, q, A))|a^a has the gauge symmetry Si^q = 
{Sxq)\x-^A- One has 

SlX. (25) 

A-5>A 

Note that we do not have to define ^^A because the equation A = A{q, q) is exactly = 0. 
We continue 

= ^(^Up,.^a)- (26) 



Sr,L = 



dLx 
dq 



hq + 



A-5>A 



dLx 
dq 



hq + 



A^A 



dLx 
d\ 



5lL 



dLx ,dLx ..dLx 
^xq + -w^(>xq + ^rrt>AA 



dq 



dq 



dX 



—Fi 

dt 1^^^ 
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This concludes the proof that L is a Lagrangian with gauge symmetry. Our result is general. 
Given any regular (i.e., non-gauge) theory and a Noether constant of motion in the canonical 
formalism, one can make this constant of motion a first class constraint and construct an 
associated Lagrangian with this gauge symmetry. 



IV. EXAMPLES IN MECHANICS 



A. Enforcing a function not being a constant of motion as a constraint 

Although we are developing the theory for implementing constants of motion as con- 
straints, let us consider an example where one implements a non-constant of motion, just 
to realize in practical terms the problems that are likely to appear. Consider the stan- 
dard Hamiltonian H{q,p) = f- — h V{q'^) {q and p are d-dimensional vectors) and try 
to implement C{q,p) = q ■ p as a, constraint. Following the above instructions we get 
P{Q:Q:^) — itT'{q ~ ^q) ^-nd A is determined as A(g, g) = m^. A substitution of both 
determinations of p and A into the extended Lagrangian yields L{q, q) — ^mqMq — V{q'^), 
where M is the matrix Mjj = 5^ — This Lagrangian is singular because the Hessian 
matrix with respect to the velocities, is (up to a multiplicative constant) identical to M, i.e. 
to a projector transverse to q. Thus L{q, q) may potentially describe a gauge theory. 

The Lagrangian momenta arc defined as p = |4 = Wlq, which indeed implies the con- 

straint qr - p ~ because Wlq = identically. The canonical Hamiltonian is just ^ + V{q'^). 
So the dynamics in phase space is given by the Dirac hamiltonian Hj:,{q, p) := H{q, p)+riq-p, 
as expected. The problem in this example is that we must require stabilization of the now 
primary constraint qr • p ~ 0. We get, as secondary constraint, ^ — q^T/' (q^) ~ 0. For a 
general potential V this gives a new condition which in its turn must be stabilized again, 
and so on. We can easily end up with incompatibility. Nothing of this kind happens if we 
choose the constraint as one of the constants of motion of the theory. 



B. Enforcing a constant of motion as a constraint 



Let us work with the same example as in the previous section, i.e. H(q, p) — 2m 
but now in R^, and with C{q,p) = e'^^q^p^. The latter is nothing but one of the conserved 
angular momenta. With this we get P^{q,q,\) — m{q^ — Xe^^^q^). Insertion of P{q,q,X) 
into the constraint C determines 



a 



(27) 



with a := {q ) + {q ) - Upon evaluation and elimination of p and A we obtain from the 
extended Lagrangian Le the new Lagrangian 



L{qq) = 2" 



a 



m 



V{q') ^ ^q 



V{q') , 



(28) 
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with the projector 



,3mn,3jA;„m„j 

^nk ^ ^nk _ t ^ V ^ (29) 

a 

It is easy to check that the projector M has := e'^-^'^g-' as the nuU vector. 

Now we work with the Lagrangian (28). The Lagrangian momenta are p — dL/dq — Mg. 
The canonical Hamiltonian becomes again ^ + ¥{0^) but there is the primary constraint 
e^i^q^p^ which is now derived from the definition of the canonical momenta and the use 
of the null vector for M. Thus the Dirac Hamiltonian is H£){q,p) := H{q,p) + rje^^^q^p^. 
Stabilization of this constraint is trivial and there are no secondary constraints in phase 
space. In agreement with this fact, one can check that the Lagrangian (28) does not yield 
constraints in tangent (i.e., configuration- velocity) space. 

One can identify the gauge transformation for L as dtq^ — e(i){g*, C}\p^p^ a-*^a = 
—e{t)e^'^^qK It is more instructive to read it by taking cylindrical coordinates z,p,9] then 
5lz — 0, 5lP — 0, 5l0 — e. In these coordinates the Lagrangian (28) is 

L = \m{z' + p2) - V{z^ + p"). (30) 

Now the gauge symmetry becomes obvious because there is no dependence on the angular 
variable in the Lagrangian. Indeed the variable B is purely gauge. The original, non-gauge, 
Lagrangian was Lng = |m(i;^ + + p^^^) — V{z^ -|- p^) , so we see that the whole procedure 
boils down to getting rid of the piece (?Q'^ . This term was invariant under rigid translations 
for the variable that is, rigid rotations around the z axis. The disappearance of this term 
makes these rotations a gauge symmetry. 

An illuminating consideration can be drawn from this example. At first sight it could 
come as a surprise that the implementation of the constraint, which requires the vanishing 
of the "angular momentum" along the z axis, allows for motions whose projection to the 
x-y plane has arbitrary dependence in the variable Q. The correct way of looking at it is the 
other way around: in promoting the constant of motion e'^^q^p^ to be a constraint, we are 
also promoting it from being a rigid symmetry generator to a gauge generator; consequently, 
the rotations around the z axis are promoted to gauge transformations. In group theoretical 
terms, the implementation of e'^^qip^ as a constraint has the consequence that a subgroup 
of the original rigid symmetry 5'0(3) gets gauged, precisely that of the rotations around the 
z axis. 



C. Relativistic massive spinless particle 

Consider the Lagrangian (spacetimc indices will be mostly supressed) Lng = |mx^ in 
Minkowski spacetime with riy^y = diag(l, — 1, ■ ■ ■ , — 1), and the rest mass m. Its associated 
Hamiltonian is H = j^P^- AH the momenta are constants of motion, so we can try to 
implement them as constraints. We then get the extended Lagrangian 

Le = pi - 7r~P^ - A(p-a), (31) 
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where in component notation \{p — a) = A^(p^ — a^), and a^u is a constant 4- vector. Elimi- 
nation of the momenta yields the Lagrangian 

La = ^m{x-Xf + Xa, (32) 

which has the gauge symmetry 6xx^ — €^{t), SxX^ — ^^{t), with e'*(T) being arbitrary 
infinitesimal functions of the evolution parameter. If we further eliminate the variables A'^, 
which have by now acquired the status of auxiliary variables, we end up with the Lagrangian 

L = ax — - — . (33) 

The last term is an irrelevant constant. Note that the EOM for L are void: every trajectory 
is a solution of the EOM. This conclusion should not be surprising because all translational 
symmetries in the Minkowski target space have been gauged, which results in making any 
trajectory acceptable as a solution of the EOM. We have simply introduced too much gauge 
freedom. 

Instead of trying to gauge the rigid translations in the target space, we could have decided 
to gauge the rigid translations along the world line, that is, the rigid translations in the 
evolution parameter. Its associated symmetry in phase space is S^x — ex, 5eP — 0, with e an 
infinitesimal constant and the generator is the constant of motion Let us fix the value of 
this constant of motion so that = m? and require this relation to become a constraint. This 
value = m? selects trajectories with unit velocity in Minkowski spacetime, x"^ = 1, but 
after enforcing this constant of motion as a constraint, a very different setting emerges, as we 
will see. For later convenience we consider the rescaled constant of motion C = 7r—(p'^ — m'^)- 
In this case, 

Le = pi - -^(/-m^), (34) 

2m 2m 

and elimination of the momenta gives 

Lx = —, + -Am, (35) 

^ 2(1 + A) 2 ' ^ ' 

which indeed has the gauge symmetry ^\x — e(r)Y^, 5a A = ^(t), obtained under the rules 
given in Section III. Addition to La of an irrelevant constant m/2 (which does not affect 
the dynamics), and a redefinition A — >■ A — 1 allows to write the modified Lagrangian (for 
which we keep the same notation) as 

-^A = + \^'^^ (36) 

with gauge transformations ^\x — e(T)|, ^a-^ = Defining as a new arbitrary function 
^ — ^1 the infinitesimal gauge transformations read 

bxx = ix, 5a A = ^(^A), (37) 
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which show directly that x is a scalar and A a scalar density under the reparamctrization 
T T — ^. The reader may rightly recognize in L\ the familiar Wheeler-Polyakov's La- 
grangian [20, 21] 

Lwp = -\{e-\T)x^{T)x^{T) + e(r)m2) , (38) 

with A = —me. The auxiliary variable e(r) is an einbein (i.e. square-root of the world-line 
metric) and r is the world-line parameter (" label time"). It can be easily checked that the 
corresponding action for Lx is invariant under finite reparameterizations of the label time, 
T !->■ t' — f{T), which, in the active view of reparameterization invariance, read 

x^ir) ^ x'^ir) = x^{r\r)) , A(r) ^ A'(r) = (^^^^) A(r^(r)) . (39) 

Here /(r) is an arbitrary monotonically increasing function of r. It is easy to check that 
the finite transformations (39) can be obtained from the infinitesimal transformations (37) 
if we set /(r) = r — ^ and successively iterate. 

The next step is to get rid of the variable A via the scheme presented in Section III D. The 
final Lagrangian L becomes L — mV^, which coincides with the usual square root world-line 
Lagrangian for relativistic particle. The corresponding action is well known to be invariant 
under reparameterizations of the label time (i.e. under the first transformation in (39)). We 
have thus succeeded in making the original theory invariant under reparamctrizations (or 
diffeomorphisms) . As a bonus wc have recovered the (on mass-shell) equivalence between 
L\YP and the square root world-line Lagrangian. 



V. THE GAUGE PRINCIPLE IN RELATIVISTIC FIELD THEORY 
A. The minimal setting 

Let us apply our results to a non-abelian field theory. For definiteness we will consider 
a A^-component complex scalar field that transforms under the fundamental representation 
of SU{N) group. The corresponding (non-gauge) Lagrangian density for the free fields is 
given by 

>Cng = r^idf.fpl-idA) - m''(t>*-(t>- (40) 
This Lagrangian has clearly SU [N] rigid symmetry 

5(t> = ie«ra</>, 5ct>* = -ie''(t>*Ta, (41) 

(Note henceforth that the action of the hermitian matrix in (l>*Ta undergoes a transposi- 
tion with respect to the action of in Ta4>) with being infinitesimal constants and Ta the 
hermitian (A^ x A^) matrices spanning a basis of the Lie algebra of SU (N), [Ta, T5] = if^f^Tc. 
To make the rigid transformation gauge we proceed along the methods outlined in Sections 
II and III. Let us first move the description in phase space. The Lagrangian definition of 
the momenta is 

TT = do4>*, n* = do4>, (42) 
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and the Hamiltonian density becomes 

n = 7V*-7T + {Vi<t>*) ■ + m^4>*-(t>. (43) 

The constants of motion which generate the rigid SU{N) symmetry are obtained as coeffi- 
cients of the infinitesimal constants in the space integration of the time component of the 
conserved current, which is computed by standard Noether methods (see, e.g. [9]). We get 

fix) = ie"[7r(a:)- Ta<t>{x) - (t>* {x)Ta ■ 7r*(x)] . (44) 

The generator G — e^-Ga ■— J (Pxj^{x) indeed generates (41) together with 

Stt = -ie^TzTa, Stt* = ie^cfy^Ta. (45) 

These transformations are in full agreement with the definition of the Lagrangian momenta 
(42). The algebra of the generators 

Ga = ij d^x [7r(x) • Tact>{x) - (t>*{x)Ta ■ 7r*{x)] , (46) 

is {Ga, Gb} — —flipc- The opposite sign in front of the structure constant fl^ is a direct 
consequence of the conventional choice [T^, T^] = ifah^c- Contact with our results from 
Section II can be established by taking cj^j, = —/„;,. 
The extended Lagrangian now takes the form 

Le = y d^xC^ = j d^X (it ■ 4) + 4)* ■ TT* - IT* ■ TZ - (V^^*) • (V^^) 

- m''ct>* ■ (f> - iX^TT ■ Ta(t> - (t>*T, ■ TT*)) . (47) 

The gauge transformations for are given by (41) and (45), but with now being an 
arbitrary infinitesimal function of time, together with the analogous of (8) 

5\\x) := 9oe"(i) - f^/{t)\^{x) = {D^e{t)Y. (48) 

Next we proceed as in Section III to construct the Lagrangian C\. We obtain, after some 
simple computations 

£a = {Do<t>nDo(P) - (V,0*) • (V,</>) - m'4>*-(P, (49) 

with 

L>o0 do(l> - i\''Ta(t>, {Do(t>y := 9o</>* + iA>*7; , (50) 

being the usual gauge covariant derivatives with the standard covariance condition S{DQ(f)) = 
ie''{x)TaDo(p. 
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B. Finishing the job 

We have succeeded with C\ in implementing gauge invariance in a restricted form. In 
fact, we have implemented it in the most minimal way, by adding as many new fields - the 
old Lagrange multipliers - as dimensions of the original rigid group we have gauged, and by 
restricting the infinitesimal parameters e°'{t) of the gauge transformation to be only functions 
of time, albeit arbitrary. On the other hand, the above implementation was so minimal that 
we have lost a big chunk of the Poincare invariance along the way. Looking at the structure 
of the term -Do</> it is clear that if Poincare transformations arc to be implemented in their 
entirety, the fields A" arc nothing else than the time components of vector fields A^, as 
dQ(f) are time components of the vector fields d^cp. Now we can in a single stroke restore 
full Poincare invariance and also let the gauge parameters to have arbitrary dependence on 
all the space-time coordinates. We just need to mimic what has been done for the time 
coordinate for all the space coordinates. In this way, gauge invariance is trivially preserved 
and we recover Poincare invariance. Then the term dicf) in the Lagrangian (49) must be 
modified to Di(f) := di<p — iA^Ta<p and similarly for di(f)*. The gauge transformations for the 
gauge fields will be the generalization of (48), namely SA^^^x) := c?^e"(a;) — f^^e^{x)A1^{x) . 
All in all we end up with the well known Lagrangian 

C = rj'"'{D^<l>*) ■ {D,(t>) - m^ct>* ■ <t> , (51) 

which is the Lagrangian for the minimal coupling of the complex scalar fields with the gauge 
field. 



C. The direct way: De Donder— Weyl formalism 

The way of finishing the job in the previous subsection leaves us with the uneasiness of 
having done it with some artifice. The problem is that the standard canonical formalism de- 
stroys the explicit Lorentz invariance and the procedure in subsection V A ends up with truly 
destroying Lorentz invariance, which then must be restored "by hand", as done in subsection 
VB. Fortunately there is a better way. De Donder- Weyl formalism [18], which preserves 
manifest Lorentz invariance in phase space, is a more suited tool to do the job. Let us go 
back to the Lagrangian (40) and define the Lorentz 4-component momenta (polymomenta) 

by 

The Hamiltonian, defined in the De Donder- Weyl formalism (DWF) through tv^ ■ d^^) -\- 
TT* ^ ■ d^cff — £ng, becomes 

T^Dw = 7r'^-7r*'^r?^, + mV-</'- (53) 

To write the extended Lagrangian we will use all four components of the SU (N) conserved 
currents, = i{7V^Ta ■ <p — 4>*Ta ■ tt*^). This is the natural way in DWF to maintain a 
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manifest Lorentz invariancc [19]. The associated multipliers A^^ are then Lorentz 4- vectors. 
De Donder-Weyl's extended Lagrangian can be then written as 

Finally, applying the methods introduced in Section III, we can successively construct La- 
grangians Cx and C By calling the latter as >Cdw we obtain 

>Cdw = r^>"'{D^(f>*)-{D,(j>) - m^(f>*-(f>, (55) 

with the covariant derivatives as defined above; D^(f) :— d^cf) — iAl^Ta(l>, etc. By finding jCbw 
we have gained a new conceptual access to gauge field theories in flat space-time. 

From here on, the rest is straightforward. One can find the curvature [D^, D^], which 
transforms under the adjoint representation of the gauge group and allows for a simple 
construction of a gauge invariant Lagrangian with kinetic terms for the Yang-Mills gauge 
fields — and a bonus of new couplings in the non abelian case. With covariant derivatives 
and curvatures at one's disposal one can analogously formulate other gauge field theories 
such as Chern-Simons gauge theory or BF gauge theory [22]. Non-local gauge invariants 
like Wilson loops or effective gluon masses [23] are also at hand. 

We have worked out the case of A^-component complex scalar field transforming under 
the SU (N) fundamental representation but we could have done the same, e.g. for the real- 
valued field multiplet in the SO{N) fundamental representation and for the spinorial case 
(e.g., for Dirac or Rarita-Schwinger fields). Note that the abelian case is recovered just as 
a particular case, as it should be. 



VI. WORLD SHEET GENERAL COVARIANCE: THE NAMBU-GOTO CLOSED 
STRING 

As another relevant example, we consider the non-gauge Lagrangian 

with the world-shcct metric hab = diag(l, —1) and the target-space (or background) metric 
f]^ll' = diag(l, — 1,...,— 1). Tis the string tension. For simplicity we will in the following 
work with natural units where T = 1. The target-space functions x^{T,a) describe the 
spacetime embedding of the world sheet. In the following we will suppress the target-space 
indices. Our aim now is to gauge the world-sheet rigid translational symmetry 

S,x = e^dax. (57) 

To prevent any conflicting issue concerning the "spatial" (a) boundary conditions we will 
deal exclusively in this section with the closed string. Following Section VC, the De Bonder- 
Weyl polymomenta are p" — — h"'''dbX, and the corresponding De Donder-Weyl 

Hamiltonian becomes 

^DW = p'^daX - £ng = ^/labpV- (58) 
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The Noether conserved current associated with the symmetry (57) is found by ordinary 
methods to be 

r = e'[p''K,p'^ - \5tp%acP'^ . (59) 

In addition to we need also to know ^^p". To compute it we resort momentarily to the 
standard canonical formalism and proceed as follow. The world sheet r-component of the 
current is 

jO = e'(^'h,p' - \5yhacP'^ , (60) 
where p^ = h}^dix = —x', so has the explicit form 

j° = 6° (^{pr - 1 [{py - {x'T]) + e\p'x') = ^ [{py + {x'r] + ^^(pV) . (6i) 

Prom this expression the transformations of p° mediated by the corresponding Noether 
charge read 

S,p' = Jda'{p'{T,a),j\T,a')} = d^{e'x' + eV) . (62) 

In deriving (62) we have allowed for e" to be an arbitrary infinitesimal world-sheet function 
to prepare the formalism for the gauge transformations we want to implement. 

By rewriting with the help of De Donder-Weyls' polymomenta we get S^p^ — di (e^p°— 
e^p^). Since in the DWF all polymomenta play the same role, we infer that the general 
transformation law for p° is 

5,p'' = 9,(ey-ey). (63) 

This should be coupled together with transformations (57) which in terms of the De Donder- 
Weyl variables read 

6,x = e'^habp'. (64) 

This last transformation also naturally follows from our definition of variations (cf. 
Eq.(4)), namely 

S,x = jda' {x{r,a),j\r,a')] = e^habp', (65) 

as it, of course, should. 

Next, in order to proceed with our program, we define the extended Lagrangian with 
Lagrange multipliers Aah- By remembering that target-space indices are suppressed we 
obtain 

= p'^daX - ^habP^'p' - iBabpV, (66) 
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where ^Bab '■= — \habAcdh'^'^ is symmetric and traceless. This shows that although we 
initially had three free Lagrange multipliers {Aab is symmetric) we end up with only two, 
because of the particular structure of the current J" and the dimensionality of the world 
sheet. 

Notice the important fact that the new EOM for Ce imply daX = {hab + Bab)p^-, and 
therefore expression (64), originated from (57) before the implementation of the Lagrange 
multipliers, needs to be reformulated to S^x — e"'{hab + Bab)p^- In turn this means (cf. 
Eq.(65)) that the conserved current needs to be reformulated. It should be noticed that 
a redefinition of currents has not been requisite in the previously discussed systems (apart 
from relativistic particle in Section IV C) because the Noether currents — coming from rigid 
(target-space) symmetries, do not change when the constraints are imposed. In contrast, here 
we deal with currents that come from rigid world-sheet symmetries and these are influenced 
when we change C^g to C^- Clearly, the same scenario occurs also for relativistic particle 
discussed in Section IV C, but there the change from 5^x — ep/m to 5f^x — ep{l + A)/m can 
be assimilated into a redeflnition of e without any extra consequences. This is not the case 
here (see our discussion later on). It is also important to observe that as defined by 
Eq.(63) is not altered because the metric tensor does not appear in expression (63) and one 
can check that the changes in the current are exactly absorbed, as regards the computation 
of SeP^-i with the redeflnition of the relation between daX and p^ , already mentioned. The 
above outlined redeflnition of the conserved current is just the flrst step in an iteration 
process, with the aim of consistency, of which we know that at every step the current will 
be quadratic in the momenta. Thus this process will result in a flnal extended Lagrangian 
of the general form 

A = p''daX-]^CabP''p\ (67) 

where Cab, which we take symmetric, contains all the information about the Lagrange mul- 
tipliers. Seen in retrospect, (66) should be interpreted as the first order expansion of Cab 
around the world sheet Minkowski metric, so that Cab = hab + Bab, with the coefficients Bab 
now taken infinitesimal. Once this observation is taken into account, we note that the trace- 
lessness condition for Bab amounts to the condition det Cab = — 1 for this Cab = hab + Bab- 
Thus det Cab = — 1 is valid at flrst order around hab- Repeated iterations of the inflnitesimal 
change hab ^ hab + Bab will be expected to preserve this condition (cf. Subsection VI B). 
Thus we end up with the result that the flnal extended Lagrangian is supplemented by the 
condition 

detC„6 = -1. (68) 

The consequences of (68) will be explored later on, in the next subsection. 

If our inputs are correct, the Lagrangian (67) should exhibit gauge freedom under the 
transformations (with arbitrary inflnitesimal functions), 

5,x = e"C„6p^ = dbieV - e'^p") , (69) 

and a certain (so far unknown) transformation 6^Cab- This means that S^^Cab should be such 
that together with (69) it should leave the Lagrangian C^i quasi-invariant, i.e,, with deCf 
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being a divergence. Let us now prove the consistency of our scheme by providing the exphcit 
form for 5^Cab- To this end we first write 

5,C, = i6y)daX + p''dai5,x) - Cab{S,p'')p'' - ^(^eC^^^V, (70) 

and notice that the first term is aheady a divergence because 

{6,p'')daX = db{eY - e''p^)daX = ^(.[(ey - ey)4a;] . (71) 
Thus ((div.) stands for divergences), 

S,Ci - (div.)+p"a„(e^Ce6/)-C„,(5,pV-^(5,C„;,)py 

= (div.) +p-da{e''Ca)p''+p''e'Ca{dap'') - C„6(9e(ey - eV))p' - ^(<^eC„6)py 

= (div.) +p''da{e'C,b)p' + p''e''C,b{daP') + (eV - 6^)9,(^6/) - ^(5eC„,)py 
= (div.) +p«9„(e^Ce5)p' +p»e'^Ce6(9„p^) + (eV - eV)(9cC„;.)p^ 

+ {e'p'' - e^p')Cab{dcp'') - l{SeCab)p''p' . (72) 

Consider the next to the last term in (72), i.e., (e"^p" — e"'p'^)Cab{dcP^)- The second piece 
cancels another term in (72), whereas the first piece can be written as 

ey^^a6(5y) = \e'^Cabdc{p''p') = (div.)-^ae(e'=C„,)py. (73) 
All in all we end up with 

= (div.) + ^p"(9„(e^Ce6) + d,{e'^cj)p' + (e^ - eV)(9cC„6)p' 
- \dc{e'Cab)P^p' -]^{5,Ca^py' 

= (div.) + \p''(da{e^C,b) + d,{e^C,,S)p' + p\'{d,Cah)p' 

-lp-(^e'^(daC,,) + e'^(d,C,a))p'' - ^5e(e^C„6)py - \{SeCa,)p''p\ (74) 
which implies that under the transformation 

S.Cab = e'^dcCab + Ccbdae' + Ca.dbe'' - Cabdc€% (75) 

the Lagrangian £f is indeed quasi-invariant. Note that this solution (75) for the transforma- 
tions rules of Cab is unique. Equation (75) is the Lie derivative of a covariant tensor density 
(0,2) of weight —1 along e, i.e., S^Cab — £eCab- Its inverse matrix, which we denote as 
C"*^ will then be a contravariant tensor density (2, 0) of weight +1, which then transforms 
according to 

5,^^ = e'^dcC^ - C^d^e" - C^d^e^ + C^d^e" = £^0"^ (76) 
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The result (75) is a very good news because the ehmination of the momenta from their own 
EOM in (66) produces the Lagrangian C\ — which in this context is more reasonable to 
denote as Cq (and similarly substitute 5\ by 5^) — which reads 

Cc{x,daX,Ci,c) = Ci{x,daX,P''{x,daX,Cde,),Cde) = ^^C^daXd^X . (77) 

The latter is a scalar density under the transformations (57) and (76), indeed 5^Cc — 
da{e"'Cc). Because transformations (57) and (76) are respectively Lie derivatives for scalars 
and for tensor densities, they — similarly as in the general relativity [9] — express diffeo- 
morphism invariance (or general covariance) of the theory. 

A. The condition det C = — 1 

The Lagrangian Cq is not the end of the story because the auxiliary variables Cab satisfy 
the additional condition det Caf, = —1. First notice that this condition is compatible with 
the gauge symmetry because det Cab behaves as a scalar under the gauge transformation 
(75), 

5ActC = e'^daidctC). (78) 

As a by-product wc see that by requiring the extra constraint det Cab = —1 the gauge 
freedom stays intact. 

In practice, one may consider two ways to implement the condition detCab = — 1 into 
jCc- One possible procedure is to introduce new gauge freedom by defining Cab — ;^/= Qab, 
with Qab an arbitrary symmetric tensor in the world sheet of signature {+,—}, and g :— 
detg^aft (note that det Qab) = — 1 and C"'' = y/~99^^)- The new gauge freedom is Weyl 
invariance, gab A(r, a)gab- This new gauge freedom compensates for the fact that gab has 
three components whereas Cab had only two. The result is the familiar non-linear a model 
Lagrangian [24-26] for bosonic string theory, 

= ^V^g'^'daxdax. (79) 

It is well known that at the classical level one can eliminate gab, which are an auxiliary 
variables in (79), by plugging their own EOM into (79). The result is the Nambu-Goto 
Lagrangian. Quantum mechanically is the issue more delicate. Instead of eliminating gab 
via its EOM, one should perform a Feynman path integral, and use the standard Fadeev- 
Popov procedure to deal with the local symmetries and gauge fixing. When this is done 
correctly [21], one finds that there is a conformal anomaly unless the target-space dimension 
is D — 26. But even in 26 dimensions it is not yet clear whether off mass-shell fluctuations 
of the Nambu-Goto and the non-liner cr-model actions contribute in the same way, say into 
string partition function. As we are interested here only in classical level description we will 
not push this point further. 

The second procedure consists in enforcing det C = —1 with a Lagrange multiplier. One 
modifies the Lagrangian (77) so that the new Lagrangian is 

JCc ^ I C'daXdbX + X{t - 1) , (80) 
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where t := \/— dct Cab (the square root is introduced for a technical convenience). Since the 
first term in (80) is already a scalar density, the transformation properties of the multiplier 
A must be also those of a scalar density, i.e. 6^X = (9a(e"A). Using the fact that Cab have 
become auxiliary variables for (80), we obtain from their own EOM that Cab = jidaxdbX, 
and therefore t is determined as 

t = -^{—det daxdbx)^ . (81) 

Plugging this result into (80) we get 

Ic ^ ^^/Xi-detdaxdbx)-* - X. (82) 

Now the multiplier A has turned an auxiliary variable. Its EOM determines A = 
{— det daX dbx)2 . Substitution of this result into (82) yields the Nambu-Goto Lagrangian 

>Cng = (- det daX dbx) ^ . (83) 

This again reconfirms the fact that on mass-shell = >Cng- 



B. Further considerations 

There is a strong parallelism between our way of obtaining the world-sheet general covari- 
ance and the approach [27] to general relativity out of the requirement of self consistency 
of the coupling of the energy momentum tensor of an initially Minkowskian theory to a 
massless spin-2 field. The presence of the coupling term produces changes in the energy 
momentum tensor which in its turn redefine the coupling term, making it non linear in the 
spin-2 field. An interaction procedure is set to work and the final result is the appearance 
of the metric tensor field and general covariance. In our case the Lagrange multipliers Bab 
play the role of the spin-2 field. A self consistency requirement also appears because the 
conserved current for world sheet translation invariance has changed due to the presence of 
the new term with the multipliers. In fact in the DWF we enforce all the components of the 
current to become constraints, and thus the Lagrange multipliers Bab are in fact coupled 
to the energy momentum tensor. The difference is that in our case, due to the particular 
structure of the current, we end up with a density tensor field Cab of weight —1 that must 
satisfy det Cac = — 1- 

Let us elaborate a bit more on the requirement det C^^ = —1. This condition is crucial 
for our purposes. In fact we have found the fulfillment of this condition for configurations of 
Cab around the flat spacetime metric and we have checked that the extension of this result 
to any configuration is fully compatible with gauge freedom. We could also argue that since 
we have found only two degrees of freedom — those of traceless symmetric Bab — around 
the flat spacetime metric, to preserve this number we must accept that the components of 
Cab are constrained by a condition of the type fiCab) = constant. If we make the reasonable 
assumption that this condition is geometrical — since the Lagrangian (83) already is — , 
we conclude that it should be a scalar under diffeomorphisms. But the only scalar we can 
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construct out of the components of the tensor density Cab is just its determinant, and to fix 
its value we need only to consider the configurations around hab- 
it is remarkable that as a way to perform the covariant quantization of the bosonic string, 
Kato and Ogawa [28] used essentially the Lagrangian (80) as a Lagrangian equivalent to (79). 
On the other hand, Siegel [29], see also [30], used the extended Lagrangian (67) with the 
specific requirement detCa^ — —1. In our approach (67) and (80) are consequences of 
gauging the world-sheet rigid translational symmetry of the Lagrangian (56). 

Finally let us stress that the dimensionality of the world sheet plays a crucial role in 
our derivation of the Nambu-Goto Lagrangian (83) through gauging the rigid world-sheet 
translational symmetry (57). It is only when the world sheet is 2-dimensional that the 
Lagrange multipliers are constrained so as to satisfy an additional condition which eventually 
leads to the requirement det Cab = ~1- 

VII. CONCLUSIONS 

Let us summarize our findings. Our starting point is a non-gauge theory, defined by a 
regular Lagrangian Lng. We assume that in the phase space formulation such a theory has 
a Lie algebra of time independent constants of motion. Next we enforce these constants of 
motion as first class constraints by adding them to the Hamiltonian with a set of Lagrange 
multipliers. Then we perform the inverse Legendre transformation to end up with a new 
(extended) Lagrangian Le whose configuration space now includes the Lagrange multipliers 
as new variables. We then observe that this new theory has gauge symmetries and that the 
gauge group is generated by the constraints, as expected. We also observe that in general the 
new variables are auxiliary and that they can be further eliminated from the formalism by 
plugging into the new Lagrangian their determination through their own equations of motion. 
This yields the final gauge invariant Lagrangian L. This last step may be problematic with 
regard to quantization because the final theory will in general be of non-polynomial nature. 
Another option is to enlarge the theory with the addition of new gauge invariant terms 
that make these auxiliary variables dynamical. The passage from L^g to L, is schematically 
illustrated in the sequence diagram in Fig.l. 

In the special case of relativistic field theories we have noticed that our program is best 
carried out if the canonical setting is taken along the lines of the De Bonder- Weyl approach. 
Such formalism is particularly suitable because it keeps manifest Lorentz invariance from 
the very scratch. The simplicity with which this gauging procedure can be performed within 
this formalism is remarkable. 

We have illustrated the DWF by applying it to the case of A'"-component complex scalar 
field transforming under the SU (N) fundamental representation but we could have done the 
same, e.g. for the real- valued field multiplet in the SO{N) fundamental representation, for 
the spinorial case, etc. It should be, nevertheless, noted that the role of the DWF is purely 
instrumental, and that once the Lagrangian for the gauge theory has been obtained (see, 
for instance, Eq.(55)), one can proceed either with Lagrangian or with standard canonical 
methods, without having to rely again on the DWF. 

As another relevant example we have derived the Nambu-Goto Lagrangian for the closed 
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^ng.(g.,i,p) =pq-H{p,q) 
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FIG. 1: The sequence diagram summarizing the basie logieal steps leading from Lng to L. The abbrevia- 
tion LM stands for Legendre map while F denotes some phase-space function which is linear in e and its 
derivatives. 

bosonic string by gauging the world sheet rigid translational symmetry of a non-gauge 
Lagrangian. Our strategy has again rehed on the DWF and it entailed an iteration procedure 
very close in spirit to the approach to Einstein's general theory of gravitation [27] in which 
a consistency argument on the coupling of a massless spin 2 field with the total energy- 
momentum tensor (including matter fields) yields ultimately the Einstein-Hilbert action. 
It should be, however, stressed that because in our reasonings the dimensionality of the 
world-sheet has played a crucial role, it is not yet clear if a similar iterative procedure can 
be applied, e.g., to relativistic Dirac-Nambu-Goto membranes (or p-branes). 

The above considered examples clearly indicate that the gauge principle, i.e., the gauging 
of a rigid group of symmetries, can be alternatively recast in the language of constrained 
systems with the gauge fields appearing first as Lagrange multipliers for the enforcement of 
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the constants of motion as constraints. The rationale of the procedure is based on the fact 
that rigid symmetries are generated by constants of motion, whereas gauge symmetries by 
first class constraints. Thus to gauge a group of rigid symmetries is tantamount to enforce 
the generating constants of motion as constraints. Note also that the role of the gauge 
fields as multipliers is temporary, because after the implementation of the inverse Legendre 
transformation they typically become auxiliary variables. Finally, when the Lagrangian is 
modified with new gauge invariant additions to provide for kinetic terms for the gauge fields, 
they become dynamical variables on their own. 

We notice also that the constraints Ca, directly originated from the former constants of 
motion of the non gauged theory, are primary constraints, but that docs not mean that our 
framework is limited to this kind of constraints and can not give rise to secondary constraints. 
On the contrary, the examples provided in section V show that, due to the presence of the 
kinetic terms for the gauge fields - which are the former Lagrangian multipliers - in the 
final Lagrangian, secondary constraints may arise, as it is indeed the case for the Yang-Mills 
gauge theories. 

With the benefit of hindsight, we observe that this route of enforcing constants of motion 
as constraints could have been taken from the very beginning as an alternative way to 
the gauge principle, because at the time when the Yang-Mills theory was formulated, the 
foundations and development of the theory of constrained systems were already in place. 

We beheve that the presented formulation can be also conveniently apphed in the 't Hooft 
program [31] where the extended Lagrangians (3) formulated with the help of constants of 
motion have played a pivotal role in construction of emergent dynamical systems [32, 33]. 
This issue would deserve further investigation. 
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